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ABSTRACT 


The  factors  involved  in  the  forming  of  point -contact  tunnel 
diodes  by  the  application  of  brief  electrical  overloads  are  studied. 
As  a  test  of  the  theory,  a  diode  fabricated  from  N-doped germanium 
and  an  aluminum  catwhisker  is  used  through  the  paper. 

After  a  brief  review  of  the  thermal  and  electrical  conductiv¬ 
ities  of  heavily  doped  germanivim  at  elevated  temperatures,  an  eval¬ 
uation  is  made  of  the  critical  breakdown  voltage  (the  applied  voltage 
below  which  no  forming  of  any  type  occurs),  assviming  constant  elec¬ 
trical  and  temperature -dependent  thermal  conductivity. 

Junction  effects  are  taken  into  account  by  considering  the 
"cold"  diode  characteristic  in  predicting  a  voltage  below  which  cui> 
rent  is  limited  by  junction  effects  and  above  which  spreading  resis¬ 
tance  dominates.  This  yields  a  value  of  series  voltage  which  must 
be  subtracted  from  the  applied  pulse  voltage  in  calculating  an  effec¬ 
tive  critical  voltage.  The  critical  voltage  for  the  sample  diode  is 
calculated  and  the  result  agrees  well  with  experiment. 


For  voltages  higher  than  the  critical  voltage,  it  is  assumed 
that  a  eutectic  region  is  fornned  under  the  point,  extending  into  the 
body  of  the  semiconductor.  The  extent  of  this  region  is  computed, 
taking  into  account  thermoelectric  effects.  On  the  basis  of  this  a 
thermal  time  constant  is  computed  which  agrees  well  with  experi¬ 
ment. 

A  heat  flow  pattern  is  computed,  and,  on  the  basis  of 
this,  two  different  types  of  diode  structures  emerge.  The  first 
is  characterized  by  a  metal- semiconductor  interface  lying  in  the 
original  plane  of  the  semiconductor;  this  configuration  does  not 
lead  to  expitaxial  regrowth  of  germanium  onto  the  original  lattice 
and  exhibits  a  typical  metal-semiconductor  I-V  characteristic. 

The  second,  which  corresponds  to  heavy  forming,  is  character¬ 
ized  by  a  eutectic  region  extending  deep  into  the  body  of  the  semi¬ 
conductor;  here  the  heat  flow  pattern  is  such  that  regrowth  of  a 
single  crystal  P-N  junction  is  to  be  expected,  a  situation  which 
leads  to  the  formation  of  tunnel  diodes. 

In  a  typical  example,  solid-state  diffusion  is  found  to  be 


negligible. 


I. 


INTRODUCTION 


i 


The  object  of  the  present  paper  is  to  investigate  the  processes 
involved  when  a  point-contact  tunnel  diode  is  "formed"  by  the 
application  of  a  short  electrical  overload. 

Forming  techniques  used  in  diodes  and  transistors  have,  for 
the  most  part,  been  derived  by  trial  and  error,  principally  because 
of  the  complicated  nature  of  the  problem.  One  of  the  most  rigor¬ 
ous  treatments  in  the  literature  is  that  of  A,  C,  Sim  (1),  Sim 
solved  the  Maxwell's  equations  and  the  heat  flow  (with  source) 
equation,  making  a  number  of  assumptions  primarily  valid  for  a 
diode  of  the  type  commonly  used  for  detectors,  harmonic  generators 

$ 

and  the  like.  He  assumed  that  the  semiconductor  material  was 
nearly  intrinsic,  that  the  point  material  was  more  refactory  than  , 

the  semiconductor  and  that  "thermal"  impurity  states  (presumably 
due  to  indiffused  copper  atoms  or  lattice  dislocations)  are  the 
primary  cause  of  the  formation  of  a  P-N  junction  in  the  neighbor¬ 
hood  of  the  contact. 

None  of  these  assumptions  are  entirely  valid  in  the  present 
case.  With  the  heavily  doped  materials  used  in  tunnel  diodes,  the 
extrinsic  conductivity  cannot  be  ignored,  at  least  at  the  lower 
temperatures,  and  the  point  materials  often  melt  at  temperatures 
low  compared  to  semiconductor  melting  points.  The  P-N  junction  ^ 

is  formed  by  indiffusion  of  dopants  contained  in  the  wire,  these 
dopants  ordinarily  being  column  III  and  V  elements  in  the  case  of  the 
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elemental  semiconductors  and  column  II-VI  elements  for  the  co¬ 
valent  semiconductors.  Sim  came  to  several  conclusions  as  a 
result  of  his  analysis,  one  of  which  was  the  desirability  of  a  high- 
impedance  pulse  supply.  Burris  and  Tramburolo  (2),  as  well  as 
the  author,  on  the  other  hand,  have  found  that  even  for  relatively 
low-current  units,  a  pulse  supply  having  an  internal  impedance  of 
an  ohm  or  less  was  usually  necessary. 

Torrey  and  Whitmer  (3)  made  important  contributions  to  the 
knowledge  of  the  fields  and  thermal  gradients  in  point  contact 
diodes  in  the  course  of  their  investigation  of  breakdown  in  detector 
diodes  due  to  T -R  tube  leakage.  Since  they  were  primarily  interested 
in  catastrophic  effects,  they  did  not  delve  into  resulting  structures, 
however. 

There  was  a  reawakening  of  interest  in  the  problem  with  the 
advent  of  the  point  contact  transistor,  due  to  observed  improve¬ 
ments  in  performance  when  the  collector  was  subjected  to  forming 
pulses.  A  number  of  papers  were  written  discussing  this  phenomena, 
most  of  the  work  being  either  experimental  results  or  theoretical 
work  of  a  qualitative  nature  (5,  6,  7),  This,  however,  was  only  a 
temporary  interest,  as  the  point-contact  transistor  declined  in 
importance  as  junction  -  transistor  performance  rapidly  improved. 

There  has  been  considerable  interest  (5,  8,  9)  shown  from 
time  to  time  in  the  influence  of  surface  effects  on  formed  point 
contacts,  a  subject  which  we  shall  neglect  on  the  assumption  that 
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such  effects  are  of  relatively  small  significance  in  a  low-impedance 
device  such  as  a  tunnel  diode. 

Most  of  the  work  in  tunnel  diodes  has  centered  about  diodes 
fabricated  using  dot-alloy  techniques  (10,  11,  12),  In  these  devices, 
the  frequency  of  cutoff  may  be  increased  by  decreasing  the  junction 
area,  since  the  spreading  resistance  increases  inversely  as  the 
contact  radius,  while  the  junction  capacity  varies  as  the  square  of 
this  radius.  While  there  are  techniques  for  reducing  the  size  of 
the  contact  area  in  alloy  junctions,  it  is  difficult  to  achieve  in  this 
manner  areas  as  small  as  those  easily  obtainable  using  a  point- 
contact  configuration.  Burris  and  Tramburolo  have  done  considerable 
work  with  the  point  contact  configuration,  and  have  reported  funda¬ 
mental  oscillations  at  103  K  Me  (2,  13,  14),  Burris  and  Tramburolo 
have  found  a  wide  variation  in  the  forming  techniques  necessary  to 
form  the  diode  properly. 

The  problem  is  an  exceedingly  complex  one,  involving  at 
many  points  effects  which  are  rather  poorly  understood  at  the 
present  time.  For  this  reason,  the  present  paper  will  not  attempt 
to  cover  the  topic  in  full  generality.  Instead,  an  attempt  will  be 
made  to  take  a  particular  metal- semiconductor  combination  and 
treat  it  theoretically  and  experimentally,  it  being  anticipated  that 
some  of  the  results  so  obtained  might  serve  at  least  as  starting 
points  in  the  analysis  of  other  combinations. 
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The  wire  chosen  is  aluminum.  The  use  of  aluminum  avoids 
the  difficulty  arising  in  the  use  of  solvents  such  as  tin  or  lead, 
namely  that  of  finding  and  interpreting  phase  diagrams  relating  to 
three  elements.  In  order  to  render  the  aluminum  relatively  easy 
to  work,  1%  of  Boron  was  added.  The  effects  of  this  small  addition 
will  be  considered  later. 

The  semiconductor  chosen  was  germanium,  since  there  is 
probably  more  tabulated  data  related  to  it  than  to  any  other  semi- 
condutor.  The  germanium  was  doped  to  3  X  10  aresnic  atoms/cc, 
it  being  difficult  to  pull  reasonably  uniform  single  crystals  at  high¬ 
er  doping  levels. 

In  this  paper,  the  experimental  and  theoretical  work  will  be 
presented  simultaneously,  since  the  numerous  approximations 
involved  require  justification  before  proceeding  to  extension  of 
the  theory  built  upon  these  approximations. 

In  order  to  carry  through  an  analysis  of  the  electrical  and 
thermal  effects  in  the  neighborhood  of  a  point  contact,  it  is 
necessary  to  have  some  notion  of  the  manner  in  which  the  electrical 
and  thermal  conductivities  of  the  heavily  doped  material  varies 
with  temperature,  especially  at  high  temperatures.  Consequently, 
Section  II  will  be  devoted  to  a  review  of  the  theoretical  and  ex¬ 
perimental  estimates  of  the  thermal  and  electrical  conductivity  of 
heavily  n-doped  germanium  in  the  region  between  room  temperature 
and  the  melting  point  of  germanium. 
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The  range  of  voltages  ordinarily  used  in  forming  the  diodes  is 
such  that  junction  effects  may  be  reasonably  expected  to  play  a  part 
in  limiting  the  current.  For  this  reason.  Part  III  will  be  devoted 
to  a  theoretical  review  and  experimental  study  of  the  forward  and 
reverse  I-V  Characteristics  of  the  diode. 

It  appears  that  there  exists  a  critical  pulse  voltage,  below 
which  the  parameters  of  the  diode  undergo  no  change;  this  obviously 
must  be  the  case  if  the  diode  is  to  function  at  all.  In  Section  IV 
this  critical  voltage  will  be  calculated  and  the  theoretical  predictions 
checked  experimentally. 

In  Section  V  the  thermal  profile  of  the  diode  for  pulse  voltages 
greater  than  the  critical  voltage  is  calculated,  Estimates  are  made 
on  the  thermal  time  constant. 

It  is  of  some  practical  importance  to  be  able  to  control  the 
current  level  of  the  finished  diode.  The  factors  involved  in  this 
are  discussed  in  Section  VI  together  with  relevant  experiment 
results. 

In  Section  VII  the  information  obtained  in  Section  V  is  used 
to  calculate  the  impurity  profile  in  the  finished  diode.  This  may 
then  be  compared  to  experimental  results. 

Section  VIII  will  be  a  review  of  the  preceding  material  as 
well  as  a  short  discussion  of  some  of  the  factors  omitted  in  the 


treatment. 
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II.  TRANSPORT  PROCESSES  IN  HEAVILY  DOPED 
SEMICONDUCTORS 

A.  Thermal  Conductivity 

In  this  section  we  shall  discuss  the  thermal  conductivity 
of  semiconductors  in  general  and  germanium  in  particular. 

Figure  1  shows  the  results  of  thermal  conductivity 
measurements  carried  out  by  Abeles  (15)  and  Slack  and  Glassbrenner 
(16).  Abeles  used  a  longitudinal  heat-flow  apparatus,  while  Slack 
and  Glassbrenner  used  a  radial  flow  configuration,  in  which  a  heat¬ 
ing  wire  ran  through  the  center  of  the  sample  and  the  thermocouples 
were  placed  radially  out  from  the  center.  This  latter  method  has 
the  advantage  that,  aside  from  negligible  end  effects,  radiated  heat 
need  not  be  taken  into  account;  for  this  reason,  it  probably  yields 
the  more  reliable  data.  Unfortunately,  this  method  was  used  only 
on  intrinsic  material.  Abeles,  however,  measured  conductivities 
for  heavily  doped  materials  (0.001  ohm-cm).  At  low  temperatures 
these  conductivities  coincide  with  his  own  measurements  for  intrinsic 
material,  while  at  high  temperatures  they  lie  between  his  intrinsic 
measurements  and  those  of  Slack  and  Glassbrenner. 

The  general  picture  for  both  doped  and  undoped  material 
-1,2  o 

seems  to  be  a  region  of  T  '  dependence  from  300  K  to  around 
900°K,  followed  by  an  abrupt  upswing. 


FIGURE  1  THERMAL  CONDUCTIYITY  OF  GERMANIUM  AS 
A  FUNCTION  OF  TEMPERATURE. 
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The  following  discussion  is  not  intended  to  be  a  tutorial 
section  on  thermal  conductivity;  such  would  be  beyond  the  scope 
of  the  paper.  Rather,  it  is  intended  to  serve  as  a  guide  to  any 
extrapolations  which  might  be  necessary  in  applying  the  rather 
meagre  experimental  data  to  other  forming  problems. 

The  most  important  concept  to  grasp  in  considering  thermal 
conduction  in  a  semiconductor  is  that  of  phonon  scattering.  If  we 
consider  a  crystal  to  be  perfectly  uniform,  and  if  we  assume  that 
the  forces  binding  the  atoms  together  vary  linearly  with  the  atomic 
displacement  from  the  lattice  site,  and  if  we  further  postulate  a 
particularly  convenient  type  of  cyclic  boundry  condition,  it  is 
possible  to  solve  the  equations  of  motion  for  the  atoms  in  the  crystal 
by  deriving  a  set  of  normal  coordinates  (17,  18,  19).  When  our 
equations  of  motion  are  written  in  terms  of  these  new  coordinates, 
they  resemble  a  set  of  equations  for  G  independent  three-dimensional 
oscillators,  where  G  is  the  number  of  unit  cells  in  the  crystal.  The 
normal  coordinates  are  identified  by  a  wave  number  k.  The 
displacements  of  the  atoms  in  the  lattice  may  be  written  in  terms  of 
k,  in  such  a  manner  that  they  resemble  a  superposition  of  traveling 
waves,  (see  Equation  (1)  ). 

U(x)  =  — ^  y  e  .  b.  (k)  exp  ik  ■  X  +  iuu.  (k)  t  (1) 
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The  vibrational  energy  of  the  crystal  is  equal  to  the  sum  c£ 
the  energies  of  the  independent  oscillators.  Quantum  mechanical 
considerations  require  that  the  energy  of  each  of  the  independent 
oscillators  must  be  an  integer  multiple  of  huu,  where  uj  is  the  angular 
frequency  of  oscillation  of  the  oscillator.  Moreover,  it  is  possible 
to  superimpose  modes  whose  wave  vectors  lie  within  a  range  Ak, 

This  will  give  rise  to  a  packet  of  spatial  extent  Ax  ~  Ak  ,  This 
localized  packet  is  called  a  phonon. 

The  phonons  carry  a  heat  current  proportional  to  their 
number  N,  energy  huu ,  and  group  velocity  ^  : 

Q  =  ^  N(k)  hu)  .  (2) 

k 

A  study  of  Equation  (2)  indicates  that,  if  the  distribution  of 

-P- 

modes  in  k  space  is  symmetrical,  there  will  be  no  heat  flow.  More¬ 
over-,  if  there  is  an  asymmetry  in  the  distribution  an  no  mechanism 
for  removing  the  asymmetry,  there  will  be  a  constant  flow  of  heat. 
Since  this  is  equivalent  to  an  infinite  conductivity,  there  must  be 
mechanisms  tending  to  redistribute  energy  among  the  phonons  in 
such  a  way  as  to  make  N  symmetric  in  k. 

Our  initial  formulation  left  no  possibility  for  the  interchange 
of  energy  among  the  phonons,  so  that,  in  order  to  arrive  at  a  finite 
conductivity,  it  is  necessary  to  consider  a  less  ideal  model.  There 
are  three  possible  approaches: 
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1.  A  reconsideration  of  the  boundary  conditions  might  be 
made,  replacing  the  mathematically  convenient  by  physically  un¬ 
realistic  cyclic  conditions  with  more  a  precise  representation  of 
the  boundaries. 

2.  A  relaxation  in  the  original  requirement  that  the  crystal 
be  perfect.  Thus,  the  presence  of  vacancies,  interstitials,  sub¬ 
stitutional  impurities,  faults,  conduction  electrons,  and  other 
deviations  from  perfect  periodicity  could  cause  a  scattering  of  a 
phonon,  with  a  consequent  variation  in  the  distribution  of  phonon 
energies . 

3.  A  consideration  of  the  possibility  that  the  restoring  forces 
acting  on  the  atom  are  not  proportional  to  displacement, 

Boundary  effects  have  been  shown  to  be  negligible  down  to 
liquid  helium  temperatures  (20). 

In  considering  other  types  of  scattering  by  stationary  defects, 
we  must  divide  the  temperature  range  into  two  parts,  the  dividing 
point  being  0^^,  the  Debye  temperature.  The  Debye  temperature  is 
derived  by  considering  a  solid  to  be  a  continuum  and  taking  its 
atomic  nature  into  account  by  limiting  the  possible  number  of  lattice 
vibrational  modes  to  three  times  the  number  of  atoms  in  the  crystal. 
If  we  call  this  highest  permitted  frequency  y^,  the  is  defined  by 
Equation  ( 3 ) , 

_ 

^D  K 


(3) 
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where  K  is  Boltzmann's  constant.  0p  is  the  point  at  which  the  lattice 
heat  capacity  begins  to  level  out  (See  Figure  2),  According  to 
Klemens  (17),  the  portion  of  the  thermal  resistance  due  to  the 
scattering  of  lattice  waves  by  static  impurities  varies  with  T  below 
the  Debye  temperature  and  is  constant  above.  This  may  be  seen  by 
referring  to  Figure  2  and  writing  the  thermal  conductivity  in  the  form 
shown  in  Equation  (4) 

K  -  ^  J  S(uj)  Mu))  dm  (4) 

where  S(m)  is  the  specific  heat  per  unit  frequency  interval.  Since 
the  net  specific  heat  is  found  to  be  integrating  S(uu )  does  not  appear 
unreasonable  that  there  should  be  a  correlation  between  S(T)  and  K(T), 

For  our  purposes,  this  dependence  of  K  on  T  below  the  De¬ 
bye  temperature  and  T°  above  means  that  it  is  possible  to  evaluate 
the  importance  of  point  impurity  scattering  by  evaluating  the  room 
temperature  thermal  conductivity  of  our  highly  doped  material  and 
comparing  it  to  that  of  instrinsic  material:  if  there  is  no  difference 
at  room  temperature,  one  may  safely  ignore  point-impurity  scatter¬ 
ing  at  higher  temperatures.  Klemens  (21)  has  also  shown  what  intuitively 
seems  clear:  that,  when  substitutional  impurities  are  randomly 
distributed  (an  assumption  ordinarily  made  with  dopants),  their 
contribution  to  thermal  resistance  varies  as  the  impurity  density. 

Since  according  to  the  resistivity  data  of  Trumbore  and  Tartaglia 
(22),  the  material  used  in  our  example  has  an  impurity  density  of 
only  two  or  three  times  that  of  the  degenerate  samples  of  Abeles(15), 
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and  since  the  data  of  Abeles  shows  no  indication  of  point  impurity 
scattering,  we  may  also  assume  that  this  may  be  neglected  in  our 
case. 

The  third  effect  mentioned,  that  of  non-linear  restoring 
forces  on  the  atoms,  is  the  most  important  one  (at  least  in  the  case 
of  germanium  and  silicon)  over  a  wide  range  of  temperatures.  The 
theory  is  rather  complicated,  and  has  been  worked  out  in  detail  by 
Klemens  (17);  in  brief  the  non -quadratic  energy  terms  arising  from 
the  nonlinear  restoring  forces  give  rise  to  interactions  between 
three  different  vibrational  modes.  This  is  called  a  three -phonon 
collision,  and  can  occur  in  two  different  ways.  If  the  collision 
conserves  crystal  momentum  hk  it  is  called  a  normal  of  N-process; 
if  it  only  conserves  momentum  to  within  a  vector  of  the  reciprocal 
lattice,  it  is  called  an  umklapp  or  U-process.  The  U-processes  are 
the  ones  which  contribute  directly  to  thermal  conductivity.  At 
temperatures  above  umklapp  processes  lead  to  a  T  dependence 
of  thermal  conductivity,  shown  in  Figure  (1)  fitted  to  room  temper¬ 
ature  thermal  conductivity.  There  does  not  seem  to  be  any  universally 

accepted  explanation  for  the  difference  between  the  experimental 

-1.2  -1 
T  '  variation  and  the  theoretical  T  variation. 

In  casting  about  for  an  explanation  for  the  increase  in 
thermal  conductivity  starting  around  940°K,  perhaps  the  first  thing 
that  might  occur  to  one  is  that  the  conduction  electrons  supplied  by 
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the  heavy  dopant  concentration  might  be  carrying  an  appreciable 
amount  of  heat.  The  contribution  of  the  conduction  electrons  to 
thermal  conductivity  is  a  familiar  problem  whose  solution,  the 
Wiedeman-Franz  law  (23),  may  be  expressed  as  Equation  (5) 

Ke  =  2a  T  (J)  .  (5) 

0  is  the  electrical  conductivity  at  temperature  T,  Tq  estimate 

the  importance  of  the  electronic  thermal  conductivity,  suppose  that 

-1  -1 

we  are  dealing  with  material  of  3000  ohm  cm  conductivity, 
which  is  about  as  degenerate  as  is  commonly  in  use  for  these  diodes. 
Then,  assuming  that  the  electrical  conductivity  does  not  increase 
as  the  temperature  is  raised,  we  get  an  electronic  thermal  con¬ 
ductivity  at  1200°K  of  5. 5  X  10  ^  watts/cm  deg,  which  is  negligible 
compared  to  the  extrapolated  value  of  the  phonon  scattering 
contributions  to  the  thermal  conductivity. 

Price  (24,  25)  has  investigated  another  carrier  effect 
influencing  the  thermal  conductivity  of  semiconductors.  The 
carriers  involved  are  intrinsic,  thermally  dissociated  electron-hole 
pairs  whose  concentration  varies  with  temperature.  Thus,  thermal 
gradients  give  rise  to  diffusion  currents,  in  which  each  hole  or 
electron  carries  with  it  an  amount  of  energy  equal  to  the  forbidden 
energy  gap  width.  Slack  and  Glassbrenner  (16)  have  computed  the 
contribution  to  the  thermal  conductivity  using  Prices  formula  and 
the  data  of  Morin  and  Miata  (26)  and  MacFarlane  et  al  (27).  Their 
prediction  for  the  thermal  conductivity  for  intrinsic  germanium  at 
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the  higher  temperatures  is  a  T  dependence  and  a  thermal  conduc¬ 
tivity  of  0.  1  watt/cm  deg,  which  is  extremely  close  to  the  extra¬ 
polated  experimental  data. 

The  question  arises  as  to  the  extension  of  these  results  to 
heavily  doped  material.  The  only  manner  in  which  heavy  doping  could 
enter  into  the  ambipolar  diffusion  process  would  be  to  decrease  the 
electron  and  hole  mobility  by  the  addition  of  an  impurity  scattering 
process.  The  degree  to  which  impurity  scattering  effects  the 
mobility  decreases  with  increasing  temperature  and  there  is  some 
data  (28)  available  which  seems  to  suggest  that,  even  at  these  high 
dopings,  the  presence  of  impurities  does  not  greatly  effect  the 
mobilities  in  the  temperature  range  900-1200°K. 

We  may  then  represent  the  conductivity  of  germanium  of  all 
doping  levels  by  Figure  3.  The  fact  that  we  may,  at  least  in  an 
approximate  treatment,  neglect  the  effects  of  doping  in  considering 
thermal  conductivity  at  the  temperatures  in  which  we  are  interested 
simplifies  the  problem  significantly.  The  actual  dopant  level  in  the 
semiconductor  is  not  uniform  during  certain  stages  of  the  forming 
process;  the  doping  is  a  function  of  the  temperature,  and  if  the 
thermal  conductivity  were  a  function  of  the  doping,  the  problem  would 
become  hopelessly  complicated. 


k,  WATTS /cm  *K 
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FIGURE  3  THERMAL  CONDUCTIVITY  OF  GERMANIUM  AS 
A  FUNCTION  OF  TEMPERATURE. 
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B.  Electrical  Conductivity 

Unfortunately  for  the  present  purpose,  there  seems  to 
be  very  little  in  the  literature  concerning  the  conductivity  of  degen¬ 
erate  semiconductors  as  a  function  of  temperature.  The  present 
section  will  discuss  the  theory  known  to  be  applicable  to  more  highly 
doped  semiconductors  and  attempt  to  justify  the  assumption  that  the 
temperature  coefficient  of  resistivity  of  our  degenerate  semiconductors 
may  be  ignored  for  the  present  purpose. 

We  shall  consider  variations  of  carrier  density  and 
mobility  separately.  The  situation  with  respect  to  carrier  density 
is  simple,  so  we  shall  consider  it  first.  "At  relatively  low  temperatures 
(between  about  300°  and  1000°K)  essentially  all  of  the  conduction 
electrons  are  supplied  by  the  ionization  of  impurity  atoms,  provided 
that  the  material  is  doped  to  degeneracy.  There  is  a  considerable 
body  of  experimental  evidence  (22,  10,  29,  32)  to  indicate  that,  at 
room  temperature,  in  semiconductors  doped  to  levels  higher  than 
10^^  carriers/cc.  ,  the  activation  energy  vanishes,  and  the  donors  are 
completely  ionized.  Hence  we  may  assume  a  1:1  relationship 
between  electrically  active  impurity  atoms  and  conduction  electrons 
up  to  the  temperatures  where  excitation  from  the  valence  bonds 
becomes  important.  It  might  be  pointed  out  that  the  term  "electrically 
active"  is  used  here  to  avoid  consideration  of  occluded  GeAs  clusters 
which  sometimes  occur  in  heavily  doped  germanium  (30). 
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At  higher  temperatures,  the  electrons  excited  up  into  the 
conduction  band  from  the  valence  band  become  important.  In  a 
rigorous  treatment,  one  would  have  to  consider  the  effect  of  the 
doping  in  order  to  calculate  the  intrinsic  electron  density,  since  the 
Fermi  level  is  determined  by  both  impurity  and  intrinsic  electron 
density.  However,  the  error  in  the  Fermi  level  resulting  from 
temporarily  ignoring  impurity  electrons  is  not  sufficient  to  greatly 
effect  the  results.  A  Semi-empirical  computation  of  the  intrinsic 
electron  distribution  was  carried  out  by  Morin  and  Maita  (26), 
working  backward  from  intrinsic  conductivity  measurements  of 
Debye  and  Conwell  (32)  and  the  theoretical  mobility  dependence 
found  by  Morin  (33).  Combining  this  with  extrinsic  carrier  con¬ 
centration  yields  a  concentration  variation  of  the  form  shown  in 
Figure  4.  For  our  purposes  it  is  clear  that  variation  of  total 
carrier  density  will  have  virtually  no  effect  on  the  resistivity. 

The  variation  of  mobility  is  a  considerably  more  complex 
subject  than  that  of  carrier  concentration  variation.  Unfortunately 
for  our  purposes,  there  seems  to  be  very  little  in  the  literature 
regarding  high-temperature  electrical  conductivity  in  degenerate 
semiconductors.  The  two  most  important  scattering  mechanisms 
are  lattice  vibrational  scattering  and  impurity  scattering. 

The  former  has  been  studied  by  Bardeen  and  Shockly  (34, 

35),  who  took  as  their  model  of  a  band  structure  that  of  the  unperturbed 
lattice,  plus  perturbations  representing  the  thermal  vibrations. 


ELECTRON  CONCENTRATION 
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After  analysing  the  resulting  scattering,  Bardeen  and  Shockley 
-1.5 

predicted  a  T  dependence  of  intrinsic  mobility  on  temperature. 

Experimental  results  (31),  however,  indicate  a  temperature 

166  Z  3 3 

dependence  of  T  for  electron  mobility  and  T  '  for  hole 

mobility.  The  reason  for  the  discrepancy  is  not  clear;  it  has  been 

suggested  (31)  that  the  many -valley  characteristic  of  E(K)  may 

explain  the  anomoly.  Figure  5  shows  the  lattice  mobility 

contribution,  plotted  against  T. 

The  basic  treatment  of  scattering  by  ionized  impurities  was 

carried  out  by  Conwell  and  Weisskopf  (36),  using  a  model  similar 

to  that  used  for  Rutherford  scattering.  Subsequent  studies  (32,  37, 

38)  have  been  mostly  refinements  of  this  work.  The  Conwell- 

Weisskopf  formula  for  the  ionized  impurity  scattering  contribution 

to  the  resistance  is  given  by  Equation  (6) 


2"/^k%^(KT)V2 
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A  plot  of  pj  versus  T  for  the  case  of  germanium  doped  to  3  X  10 
carrier s/c.  c.  is  given  in  Figure  5.  Also  given  is  a  composite 
mobility  using  a  combination  formula  of  Johnson  and  Lark-Horovitz 


(39). 


The  Conwell -Weisskopf  formula  is  essentially  a  classically 


derived  expression.  Brooks  (41)  attacked  the  problem  from  a 
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quantum -mechanical  point  of  view  and  concluded  that  Equation  (6) 

2 

m  0 

should  have  a  multiplication  factor  of  — t-tsj - r  in  the  argument  of 

N  '  Tir.K 

the  logarithmic  term.  Although  the  meaning  of  K  is  not  clear,  which 
precludes  numerical  correction  of  the  Conwell-Weisskopf  expression, 
it  is  clear  that  by  suppressing  the  denominator,  the  correction  tends 
to  prevent  the  decrease  of  p  with  increasing  temperature.  This  is 
what  we  expect,  since,  by  using  the  Fermi  statistics  instead  of 
Boltzmann  statistics,  one  considers  only  those  electrons  having 
energies  around  the  Fermi  energy,  which  in  germanium  is  only 
weaJcly  dependent  on  temperature.  Hence  thermally  induced  varia¬ 
tions  in  energy  in  the  "tails"  of  the  distribution  have  relatively  little 
effect  on  the  mobility.  Now,  evenwith  the  pessimistic  Conwell-Weiss¬ 
kopf  prediction,  the  mobility  and  hence  the  conductivity  varies  only 
a  factor  of  two  between  room  temperature  and  the  melting  point.  In 
view  of  our  expectation  that  the  variation  is  even  less,  we  shall 
make  the  assumption  that  the  resistivity  is  a  constant. 

This  assumption  is  crucial  to  the  analysis.  Without  it,  it 
would  be  necessary  to  solve  the  thermal  and  electrical  equations 
simultaneously,  each  equation  involving  nonlinear  functions  of  the 
solution  of  the  other.  By  assuming  that  the  temperature  coefficient 
of  resistivity  may  be  justifiably  ignored,  one  may  solve  first  the 
Maxwell  equation  for  potential  or  current  distribution,  then  the 
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thermal  conductivity  equation,  using  the  solution  to  the  former 
as  sources  in  the  latter. 

C.  Thermoelectric  Effects 

Although  they  are  classified  as  transport  processes, 
we  shall  defer  the  discussion  of  thermoelectric  effects  until 


Section  V, 
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III.  JUNCTION  EFFECTS 

The  purpose  of  this  section  is  to  derive  a  model  of  the  diode 
which  will  permit  computations  which  take  into  account  junction 
effects . 

Two  points  should  be  made  here.  In  all  our  studies  we 
shall  assume  that  the  point  is  given  a  very  light  current  pulse  to 
anchor  it  to  the  semiconductor  before  the  main  forming  pulse  is 
applied:  otherwise,  there  would  have  been  too  much  mechanical 
instability,  and  surface  effects  might  have  played  an  excessive 
part  in  the  starting  conditions.  Secondly,  most  of  our  experiments 
were  conducted  with  cleaved  surfaces,  with  surfaces  ground  and 
etched  with  CP^,  and  with  surfaces  ground  and  etched  with  hydrogen 
peroxide  etch.  The  principal  difference  in  results  was  that  cleaved 
surfaces  tended  to  permit  more  sensitive  determination  of  contact 
as  the  wire  was  brought  in  to  touch  the  semiconductor  than  was 
possible  with  etched  surfaces.  For  a  given  wire,  contact  pressure, 
and  forming  system,  the  three  types  of  surface  preparation  had  little 
influence  in  the  final  diode. 

Since  all  our  successful  tunnel  diodes  were  made  using 
forward  pulses,  we  shall  consider  the  forward  characteristic  first. 
By  and  large,  we  may  consider  the  forward  characteristic  as  being 
divided  into  three  distinct  regions;  the  first  dominated  by  the 
tunneling  current,  if  any,  the  second  by  the  normal  diode  current, 
and  the  third  by  the  effects  of  spreading  resistance. 


-32- 


Although  the  tunneling  current  plays  no  direct  part  in  the 
forming  process,  we  shall  discuss  it  here  in  some  detail,  in  order 
to  derive  some  criterian  for  the  impurity  profile  necessary  to 
form  a  tunnel  diode. 

In  the  first  place  it  must  be  realized  that  we  are  dealing  here 
with  a  P-N  junction  formed  by  the  alloying  of  the  point  materials  and 
not  with  a  metal- semiconductor  junction.  This  fact  has  been 
established  by  thermoelectric  probing  experiments,  (42)  and  by  the 
fact  that  alloys  suitable  for  dot-alloy  metals  also  usually  work  well 
for  point -contact  wires. 

The  requirements  for  tunneling  of  the  sort  considered  here 
are  twofold:  first,  material  on  either  side  of  the  junction  must  be 
doped  to  degeneracy,  and  secondly,  the  transition  region  must  be 
very  narrow.  To  sharpen  and  give  some  numerical  significance 
to  these  terms  a  very  brief  review  of  tunnel  diode  action  will  be 
given  here.  The  theory  of  tunneling,  particularly  in  germanium,  is 
extremely  complex  and  is  still  not  completely  understood,  (43,  44, 

45)  so  that  our  review  will  be  limited  to  an  attempt  to  obtain  a 
criteria  for  the  impurity  profile  necessary  for  noticeable  tunneling 
to  occur. 

Figure  6  shows  the  conventional  representation  for  a  tunnel 
diode  energy  level  configuration.  The  material  on  both  sides  is  doped  to 
degeneracy;  that  is,  theFermilevellies  within  the  conduction  band  on 
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the  n-side  and  within  the  valence  band  on  the  p-side.  Since  with  no 
external  voltages  applied  the  Fermi  levels  on  both  sides  line  up, 
empty  states  in  the  valence  band  sit  opposite  filled  states  in  the 
conduction  band.  When  forward  bias  is  applied,  the  bands  uncross 
until  finally  there  are  no  longer  valence  band  states  opposite 
conduction  band  states.  It  should  be  noted  that  we  are  ignoring 
band  edge  smearing  and  impurity  -  band  effects  (46),  which  become 
important  principally  in  discussing  the  valley  current. 

There  is  a  finite  probability  that  an  electron  at  the  band 
edge  of,  say,  the  P-side  valence  band  may  cross  the  so-called 
forbidden  band  and  enter  the  N-type  conduction  band  with  no  loss 
of  energy.  At  zero  bias  the  density  of  states  function  on  each  side 
of  the  junction  is  such  that  the  flow  of  tunneling  electrons  in  each 
direction  is  equal  and  no  net  current  flows.  With  a  applied  bias, 
however,  the  balance  is  destroyed  and  there  is  a  net  current  flow. 

The  tunneling  current  density  depends  on  three  factors:  the 
distribution  of  available  and  filled  levels  in  the  coincident  bands,  the 
rate  at  which  electrons  strike  the  band  edge,  and  the  probability 
that  an  electron  which  strikes  the  barrier  will  tunnel  through. 

The  distribution  of  available  and  filled  levels  is  found  by 
multiplying  the  density-of- states  function  by  the  Fermi  distribution 
function.  The  tunneling  current  due  to  a  particular  energy  is  pro¬ 
portional  to  the  product  of  the  distributions  for  both  sides. 
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The  second  factor  in  the  tunneling  current  is  the  number  of 
times  that  the  electrons  strikes  the  band  edge  under  the  influence 
of  the  built-in  field  in  the  junction.  This  picture,  as  developed  by 
Franz  (47)  and  Spenke  (48),  pictures  the  electron  as  being  acted 
upon  by  a  force  e  e  due  to  the  distortion  of  the  bands.  A  quasi - 
classical  force -acceleration  law  is  assumed  to  hold: 

F  =  =  ee  =  hk  .  (6) 

The  variation  in  k  implies  a  variation  in  the  energy  E; 
each  time  k  changes  by  K  (a  vector  in  the  reciprocal  lattice)  the 
energy  of  the  electron  goes  through  all  the  values  of  the  band.  Thus 
we  can  consider  a  frequency  of  striking  the  band  edge  as  being 
described  by  Equation  (7) 


l/t. 
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where  we  have  used  a  simple  cubic  lattice  characterized  by  a  lattice 
constant  a.  e  is  the  field  strength  in  the  junction. 

The  calculation  for  the  tunneling  probability  has  been  carried 
out  by  several  workers  (43,  44,  45,  47)  and  is  still  the  subject  of 
some  controversy.  We  shall  use  the  result  of  Kane's  (43)  treatment. 


which  yields  as  the  probability  T  per  incident  electron: 
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where  B  is  a  function  of  k  but  not  of  0 .  If  we  combine  8,  7,  and  the 
distribution  term  we  get  a  representation  for  the  tunneling  current 


density: 


J  =  C, 
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D(E)  represents  here  the  density  of  state -probability  of 
occupation  terms.  This  term  is  extremely  difficult  to  evaluate  in 
the  case  of  germanium  because  of  the  many- valley  characteristic 
of  E(k)  for  germanium  and  also  because  there  are  two  valence  bands 


in  germanium.  Kane  makes  an  admittedly  crude  estimate  of 
tn^ir^m*Er 

- 2 — 2 — ^  value  of  D  at  the  current  peak,  with  being 

12h  e  0 

the  distance  between  the  Fermi  level  on  the  n  side  and  the  band 


edge. 

We  must  now  solve  for  0.  If  we  refer  to  Figure  6,  0  is 
usually  assumed  to  be  the  slope  of  the  band  edges  in  the  middle  of 
the  transition  region.  The  distance  W  from  to  x^  is  found  by 
solving  Poisson's  equation 

V"(x)  =  -  —  C.(x)  (10) 

€ 


where  ^(x)  is  the  net  charge  density  and  0  is  the  dielectric  constant. 
The  boundary  condition  is  that  the  difference  in  the  electrostatic 
potential  on  the  n-side  and  on  the  ^  -side  must  equal  E^  n  p 
this  is  simply  the  condition  that  the  Fermi  levels  on  the  two  sides 
must  be  aligned.  Simultaneously  with  (10)  it  is  necessary  to  satisfy 
the  neutrality  condition  obtained  by  equating  the  number  of  ionized 
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impurities  to  the  number  of  filled  conduction  band  states  on  the 
n  side  and  to  the  number  of  empty  valence  band  states  on  the  p-side. 
The  solution  for  the  potential  distribution  is  not  simple  in  an 
arbitrary  doping  profile,  but  it  has  been  calculated  (43)  for  the  case 
of  an  abrupt  junction.  The  result  is  given  in  Equation  (11): 


where  n  and  p  are  the  dopant  concentrations  on  the  n  and  p  sides, 
respectively,  e  is  the  dielectric  constant,  and  V  is  the  potential 

19 

difference  between  the  two  sides.  If  we  assume  that  n=p=3X10 

carriers  per  c.  c.  and  that  V  «  ^  ^  »  Equation  (11)  yields 

_7 

a  value  of  W  equal  to  9,1  X  10  cm.  Using  this,  we  may  calculate 
the  value  of  e  that  appears  in  Equations  (6,  7,  8,  and  9).  We  ignore 
the  curvature  in  the  potential  across  the  junction  and  assume  that 
any  voltage  swing  necessary  to  obtain  the  peak  current  may  be 
neglected  in  comparison  to  E^.  We  then  write  (9)  in  the  form 

T  =  ^  r 

max  »l/2  L 
(E^) 

The  constant  is  difficult  to  calculate  theoretically,  due  to 
the  complexities  of  the  band  structure  of  germanium.  However, 
capacitance  studies  (50)  tend  to  indicate  the  the  form  of  (12)  is  valid. 
These  measurements  also  indicate  that,  for  transition  regions  as 
thin  as  we  have  considered  and  for  doping  levels  comparable  to 


*l/2_l/2 
-vm  '  W 
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those  in  our  example,  the  current  density  is  on  the  order  of  10 

z 

amperes/cm  .  Burris,  using  point-contact  diodes  and  estimating 
the  junction  area  by  microscopically  examining  the  surface  after 

5 

the  point  had  been  pulled  away,  estimated  current  densities  of  10 

amp/cm  .  The  author,  using  a  sectioning  technique  which  will  be 

described  later,  got  values  up  to  110,  000  amperes/cm^.  If  we 

take  the  highest  of  these  values  and  assume  that  it  corresponds 

-7 

to  the  minimum  thickness  of  9- 1  ^  10  cm,  we  can  obtain  the 
peak  tunneling  current  as  a  function  of  W.  At  W  =  ZOO  A.,  the 
tunneling  current  is  down  to  10  amps/cm  ,  and  may  be  completely 
ignored.  It  is  obvious  that  even  if  we  make  generous  allowances 
for  the  crudity  of  some  of  the  approximations  involved,  junctions 
substantially  wider  than  a  few  hundred  angstroms  are  unlikely  to 
lead  to  tunneling  currents. 

Let  us  now  consider  the  diode  equation  current  and  its 


influence  on  the  i-v,  characteristic.  The  diode  equation  is 


Jr,,  j  =  C,  (e  '"'^-1) 
Diode  4  '  ' 


where  and  are  constants  and  V  is  the  applied  voltage.  Theory 
predicts  (51)  that  is  a  function  of  mean  free  paths  and  diffusion  lengths 
and  that  Cg  equals  e/kT ;  the  former  requires  data  that  is  difficult  to 
obtain  and  the  latter  is  only  infrequently  valid.  We  shall  not  discuss 
the  various  theories  explaining  the  discrepancies,  but  shall  instead 
resort  to  experiment.  At  this  point  we  should  mention  the  experimental 
setup  used  in  this  study. 
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Figure  7  shows  the  jig  in  which  the  wire  and  semiconductor 
were  held  during  the  electrical  processing,  A  spring-loaded  micro¬ 
meter  adjustment  with  an  oversize  thimble  permitted  accurately 
controlled  motions  with  resolution  better  than  0,  0001",  One 
problem  which  arose  with  the  aluminum  wire  was  that  electrolytic 
pointing  or  even  well -controlled  mechanical  pointing  was  impossible, 
since  the  rapid  rate  of  oxidation  of  freshly  cut  aluminum  made  it 
necessary  to  form  the  junction  almost  immediately  after  snipping. 
Hence,  in  some  of  the  diodes  formed,  the  wafer  and  wire  were 
protected  from  the  atmosphere  by  a  Incite  enclosure,  so  that 
operation  in  an  inert  atmosphere  was  possible.  To  expose  fresh 
aluminum  a  vane  of  very  pure  aluminum  was  touched  to  the  end  of 
the  wire  and  a  voltage  of  about  ten  volts  applied  to  melt  the  wire  in 
half,  forming  a  Ball  of  very  clean  aluminum  on  the  end  (see  Figure 
8),  The  vane  was  then  swung  out  of  the  way  and  the  wire  brought 
into  contact  with  the  wafer. 

Figure  9  shows  the  system  in  its  entirety,  A  sweep  capable 

of  showing  the  forward  and  reverse  characteristics  of  devices  having 

-1  -7 

peak  current  levels  ranging  from  10  to  10  amperes  was  connected 
to  the  holder  through  a  mercury  switch  (to  minimize  transients  due 
to  contact  potentials).  Also  connected  through  a  mercury  switch 
was  a  pulse  supply,  timed  by  the  gate  pulse  from  a  Tektronix 
oscilloscope  time  base  generator.  The  circuit  of  the  pulse  amplifier 
is  shown  in  Figure  10;  it  consists  of  a  vacuum  tube  cathode  follower 
driving  a  transistor  inverter,  which  in  turn  drives  a  cascaded 
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FIGURE  8  ALUMINUM  WIRE  AFTER  FLASHING. 


-42- 


FIGURE  9  FORMING  AND  MONITORING  SYSTEM. 
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series  of  emitter  followers.  The  last  two  stages  were  2N1907 
high-frequency  diodes  in  emitter  follower  configurations  with  heavy 
feedback  (about  0,7  ohms)  between  emitter  and  input.  The  circuit 
was  capable  of  putting  out  pulses  of  5-1000  psec  duration  and  with 
pulse  height  variable  (by  means  of  a  clipper  circuit  in  one  of  the 
pre-driver  emitter  followers)  from  0  to  3.8  volts.  The  rise  and 
fall  times  were  less  than  0.  5  psec  and  the  internal  impedance 
was  0.  1  ohm.  The  circuit  was  capable  of  putting  out  15  ampere 
pulses  at  two  volts,  and  was,  moreover,  capable  of  tolerating 
microsecond  alterations  in  the  load  impedance,  which  was  a  some¬ 
what  unusual  feature  of  this  application. 

In  measuring  the  diode  characteristic,  a  freshly  snipped  or 
flashed  wire  was  brought  down  onto  the  semiconductor  until  the 
sweep  display  at  its  most  sensitive  setting  indicated  current  flow. 

It  was  then  subjected  to  some  particular  pulsing  format  and  the 
picture  of  its  forward  and  reverse  characteristics  taken  with  a 
scope  camera.  A  split  brass  form  was  then  carefully  clamped 
in  place  and  epoxy  poured  around  the  diode.  The  potted  diode 
was  then  carefully  ground  and  polished  until  the  junction  was 
exposed,  after  which  it  was  examined  and  photographed  with  a 
metallographic  microscope.  A  typical  section  is  shown  in 
Figure  11  ,  The  technique  is  hampered  by  the  variation  in  hardness 
between  the  germanium,  the  aluminum  wire,  and  the  epoxy;  abrasives 
suitable  for  giving  a  good  reflective  finish  to  germanium  tend  to  scoop 
and  undercut  the  aluminum  and  the  germanium  epoxy  interface. 
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Moreover,  the  welded  point  tends  to  set  up  a  stressed  condition 
in  the  germanium  in  the  immediate  neighborhood  of  the  point 
which  the  grinding  appears  to  aggravate,  with  the  result  that 
chipping  and  "cribbing"  often  results.  However,  usable  data 
was  obtained.  Some  measurements  were  also  made  by  pulling  the 
wire  off  and  examining  the  spot  where  it  had  been  (Figure  12),  This 
is  for  a  relatively  high  current  unit.  Note  that  there  appears  to  be  a 
relatively  symmetrical  shear  area,  probably  coinciding  with  faults 
introduced  during  forming.  The  rough  region  presumably 
corresponds  to  the  P-N  junction  area. 

Part  of  the  results  are  shown  in  Figure  13,  Diodes  3-E4, 

3-25,  are  3-26  were  formed  by  pulsing  with  relatively  high -impedance 
(100  ohm)  sources;  if  we  refer  to  Equation  (13),  these  have  a  value 
of  Cg  equal  to  7.  25.  Diodes  3-27,  3-28,  and  3-27  were  formed  by 
pulsing  with  a  low-impedance  (0,  1  ohm)  supply;  these  yielded  a  value 
of  Cg  equal  to  15.  75.  The  reason  for  the  variation  of  in  diodes  is 
rather  poorly  understood  (3). 

Now  let  us  use  some  of  our  experimental  results  to  estimate 
the  degree  of  importance  of  junction  effects.  The  incremental 
conductance  of  a  diode  may  be  found  by  differentiating  Equation  (13): 

2  ^5^ 

Cp  =  ira^  C4  Cg  e  (14) 

where  a  is  the  radius  of  the  contact.  The  spreading  conductance 
is  (3) 


4a 


(15) 
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FIGURE  12  VERTICAL  VIEW  OF  JUNCTION 
AFTER  WIRE  IS  REMOVED. 
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FIGURE  13  DIODE  CHARACTERISTICS. 
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where  p  is  the  resistivity  of  the  material.  The  forward  voltage  at 
which  the  spreading  resistance  comes  into  play  may  be  found  by 
setting  (14)  equal  to  (15)  and  solving  for  x: 


*°8i0  [ 


10  LirapC^C^ 


,-3 


The  largest  value  of  a  encountered  in  our  work  was  a  =  10  cm, 

_3 

and  the  material  was  10  ohm- cm.  An  average  value  of 

obtained  by  microscopic  examinations  of  the  type  described 
earlier  was  10  amps/cm  .  This  leads  to  a  value  of  V  equal  to 
0.  568  volts.  Figure  13  shows  the  I-V  characteristic  at  high 

.3 

current  levels  for  a  diode  for  which  a  =  10  cm.  Note  the 

sudden  break  between  0.  5  and  0. 6  volts,  corresponding  to  this 

transition.  The  maximum  voltage  at  which  junction  effects  are 

important  isV=E  +4.  since  at  higher  forward  bias  the 

g  n  p 

junction  has  effectively  disappeared.  This  will  be  around  0.  7 
volts  for  germanium.  Hence  we  may  expect  that  junction  effects 
may  be  taken  into  account  by  postulating  a  series  voltage  source 
of  between  0.  57  and  0.  7  volts,  depending  upon  the  area  of  the 
contact. 


Little  success  was  obtained  in  forming  using  reverse 
pulses,  and  the  reverse  characteristics  varied  too  widely  to  permit 
characterization  in  any  simple  and  practical  manner  .  Close 
empirical  fit  was  often  possible  to  the  Zener-like  expression  in  (18) 
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The  variation  of  and  was  too  gross  to  allow  a  very  useful 
generalization  to  be  made;  roughly,  junction  effects  appeared  to 
dominate  spreading  resistance  out  to  about  1  volt. 
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IV.  COMPUTATION  OF  CRITICAL  BREAKDOWN  VOLTAGE 

When  a  metal  point  is  applied  to  a  semiconductor  and  a 
voltage  applied  across  the  combination,  the  flow  of  current  through 
the  resistive  bulk  semiconductor  material  will  generate  a  certain 
amount  of  heat,  giving  rise  to  a  temperature  gradient  within  the 
material.  Intuition  indicates  and  more  rigorous  analysis  proves 
that  the  highest  temperature  will  be  associated  with  the  material 
immediately  adjacent  to  the  point.  If  the  applied  voltage  is  high 
enough,  either  the  wire  or  the  semiconductor  will  melt  in  the 
vicinity  of  the  interface.  From  the  point  of  view  of  forming 
studies,  the  voltage  at  which  this  occurs  is  a  critical  one,  and  it  is 
the  purpose  of  this  section  to  solve  for  it  and  to  check  the  analytical 
results  by  experiment. 

Torrey  and  Whitmer  (3)  attacked  this  problem  for  the  case 
in  which  the  thermal  and  electrical  resistances  are  constant  with 
temperature,  and  Sim  solved  it  assuming  the  electrical  conductivity 
to  vary  with  temperature  and  the  thermal  conductivity  to  be  constant 
with  temperature.  Following  the  treatment  in  Section  II,  we  shall 
assume  that  the  electrical  conductivity  is  constant  while  the  thermal 
conductivity  varies  in  the  manner  shown  in  Figure  3. 

Figure  14  shows  the  metal- semiconductor  contact.  We 
assume  that  »  a^.  so  that  the  interface  will  be  an  equipotential 
plane.  Also,  if  the  metal  is  molten  at  the  interface,  we  may  assume 
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FIGURE  14  THE  METAL-SEMICONDUCTOR  CONTACT. 
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that  the  metal -semiconductor  interface  is  an  isotherm  at  the  melting 
point  of  the  contact  material.  It  is  obvious  also  that  there  is  going 
to  be  a  high  degree  of  cylindrical  symmetry  to  the  problem.  The 
form  best  suited  to  expressing  these  conditions  is  the  oblate 
spheroidal  coordinate  system.  (Figure  15) 

In  this  system,  the  orthogonal  surfaces  may  be  written  in  the 
forms  shown  in  Equations  (18)and(19). 


2,  2, 

at.  a  (1+t  ) 


=  1 


(18) 


where 


2  2 
P  .  _ 

2, .  > .2  2t2 

a  (1-1)  a  I 


2  2,2 

p  =  X  +  y 


=  1 


(19) 


(20) 


In  addition,  there  are  planes  of  constant  azimuthal  angle  (j). 


The  matrics  are 


h^  =  a[(l  +c  ^)(l-e^)] 


1H‘ 


(21) 


1/2 


In  a  conducting  medium  with  negligible  charge  accumulation, 
V  '  T  =  0,  where  J  =  the  current  density.  Thus, 


J=  -a3VV 


(22) 
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so 


V-{-OgVV)=0  (23) 

We  how  make  use  of  the  assumption  of  constant  electrical 
conductivity  which  we  developed  in  Section  II.  This  permits  us  to 
take  the  term  in  (23)  over  to  the  left  of  the  gradient  operator, 
giving  us  LaPlaces'  equation: 

-a  V^V  =  0  (24) 

or,  in  our  particular  coordinate  system: 


8 


I?  ] 


^4^ 


8^V 


(1H^)(1-  '')2f 


=  0 


(25) 


8/8  (j)  =  0  from  symmetry,  so  that  (25)  becomes: 


8 

(26) 

Assume: 

V  = 

A(e)  B(4) 

• 

(27) 

This  permits  us 

1  to  separate  the  equation 

d 

(1-e^) 

+  n(n+l)  A  =  0 

(28) 

d 

■(1H^)^‘ 

]  -  n(n+l)  B  =  0 

(29) 
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We  shall  assume  that  the  electrical  conductivity  of  the  point  material 

is  so  much  higher  than  that  of  the  semiconductor  that  the  metal - 

dA 

semiconductor  interface  is  an  equipotential.  Hence  must 
vanish  and  so  must  n.  Then  (29)  becomes 


4;  =  0  (30) 

solving  for  B, 

B  =  C  tan"^  ?,  +  D  =  V(;)  .  (31) 

To  solve  for  C  and  D  we  set  the  potential  V  equal  to  zero 
at  C  =  -  00  and  equal  to  V^,  the  effective  applied  voltage,  at  ^  =  0. 
Then: 


Let  P 
loss. 


V  =  V  [l  -  -tan"^  cl 
OL  IT  J 


(32) 


equal  the  power  input  per  unit  volvime  due  to  ohmic  power 
Then; 


_  ,„.,,2,  r  1  svi 

P=a^(VV  )  =  o,  Le:  wJ 


4V^a 
o  s 


2  2,{,2  _  2., .  „  2. 

IT  a  (I  +C  )(1+C  ) 


(33) 


This  is  the  same  as  the  Torrey-Whitmer  result  for  power 
distribution.  We  shall  now  solve  for  the  temperature  distribution. 
We  can  write  the  heat  flow  equation,  setting  the  heat  flux  out  of  the 
material  equal  to  the  ohmic  power  input. 


V  •  KVT 


4V 

Tr^a^(C^+C^)(l+C^) 


(34) 


We  have  already  assumed  that  the  metal- semiconductor 
junction  is  an  equipotential  surface;  we  may  also  see  that  it  is  an 
isothermal  surface,  at  a  temperature  equal  to  the  fusion  temperature 
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of  the  point  material.  It  has  been  shown  (1)  that,  with  configurations 
of  this  sort,  if  the  equipotential  and  isothermal  surfaces  coincide  at 
one  surface,  they  coincide  throughout.  Hence  we  are  justified  in 
setting 


8T 


8T 


=  0 


(35) 


This  reduces  Equation  (34)  to: 

h,.h 


1  d  r  dTi_^^o'^8  i 


which  reduces  to 


Integrating  once. 


4V^a 
o  s 

~2  Z 

ir  (1+^  ) 


7  AT  4V^a  cos”^^ 

K(1H^)^  =  - 


or 


dT 


4V^a 

O  8 

■/k 


IT 


cot 


1+^ 


(37) 


(38) 


(39) 


We  cannot  integrate  this  exactly,  because  of  the  irregular 
variation  of  K  withT  and  hence  with  C  •  However,  its  numerical 
solution  poses  no  particular  problem;  this  has  been  carried  out  and 
the  results  plotted  in  Figure  16. 

We  shall  now  simultaneously  check  Figure  16  and  the 


conclusions  of  Section  III  by  a  simple  experiment.  To  do  this  we 
find  the  melting  temperature  of  the  point  material  (or  of  germaniiam 
if  it  is  the  less  refractory)  and  use  Figure  16  to  find  the  effective 
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voltage.  To  this  must  be  added  the  series  voltage  due  to  junction 
effects,  as  calculated  in  the  preceeding  section. 


z 

FIGURE  15  THE  OBLATE  SPHEROIDAL  COORDINATE  SYSTEM. 


In  our  case  the  point  material  melts  before  the  semi¬ 
conductor,  at  a  temperature  which  we  may  estimate  from  the 
aluminum -boron  phase  diagram  (53)  as  being  approximately 
920°  K.  Referring  to  the  chart  of  Figure  16,  this  corresponds  to 
an  effective  applied  voltage  of  0.  62  V  necessary  to  cause  break¬ 
down.  The  diode  used  in  Section  HI  as  an  example  was  also  used 
here.  Referring  back  to  the  calculation  of  Section  HI,  we  see 
that  the  effective  series  voltage  due  to  junction  effects  is  0.  57 
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volts.  Thus  the  applied  voltage  necessary  for  breakdown  should 
be  0.  62  volts  plus  0.  57  volts,  a  total  of  1.19  volts. 

To  measure  the  breakdown  we  apply  sweep  voltage  of 
increasingly  higher  amplitude  to  the  diode,  from  a  moderately 
high-impedance  supply.  After  each  time  that  the  diode  is  subjected 
to  a  sweep  voltage  of  particular  peak  amplitude,  the  circuit  is 
switched  to  a  more  sensitive  setting  and  the  current  at  a  forward 
bias  of  0.  1  volt  measured.  As  we  shall  see  later,  the  thermal  time 
constant  of  the  diode  configuration  is  on  the  order  of  microseconds, 
so  that  the  fact  that  we  are  using  a  60-cycle  sweep  rather  than  dc 
will  have  no  effect.  The  results  of  the  measurement  are  shown  in 
Figure  17.  The  reason  for  the  slight  descending  characteristic  at 
the  lower  values  of  sweep  voltage  is  not  known.  However,  it  is 
clear  that  a  sharp  breakdown,  prestimably  due  to  a  melting  of  the 
point,  took  place  at  1.17  volts.  This  represents  an  error  of  less 
than  2  per  cent  in  our  original  estimate  of  1.  19  volts.  The  reason 
that  the  current  after  breakdown  was  not  higher  was  that,  as  soon 
as  the  diode  began  breaking  down  and  more  current  flowed,  the 
sweep  voltage  dropped  due  to  the  series  impedance  of  the  pulse 
supply.  Thus  it  was  not  necessary  to  take  into  account  the  effects 
of  exceeding  the  critical  voltage,  a  subject  which  will  be 
considered  in  the  next  section. 


1.0 


1.08 


1.12 


1.20 


1.04  1.08  1.12  1.16 

MAXIMUM  VALUE  OF  SWEEP  VOLTAGE.  VOLTS 


FIGURE  17;  1^^^  VERSUS  MAXIMUM  FORWARD  SWEEP 
VOLTAGE. 
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V.  THE  THERMAL  PROFILE  AND  DIODE  CONFIGURATION 
FOR  VOLTAGES  HIGHER  THAN  CRITICAL 
We  saw  in  the  previous  section  that,  as  the  voltage  applied 
to  a  lightly  formed  diode  was  increased,  a  value  of  voltage  was 
finally  reached  at  which  melting  of  the  point  began.  We  shall  now 
analyse  the  case  where  the  voltage  is  higher  than  this  critical 
value . 

Figure  18  shows  a  diode  being  subjected  to  heavy  forming. 
The  heating  causes  the  formation  of  a-molten  euteCtic  mixture  of 
point  material  and  semiconductor,  the  solid  germanium  having  been 
partially  melted  out  in  a  crater  beneath  the  point. 


FIGURE  18  DIODE  SUBJECTED  TO  HEAVY  FORMING. 
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It  is  quite  difficult  to  make  any  very  intelligent  guesses 
about  the  conductivity  of  the  molten  eutectic,  there  being  relatively 
little  data  available  on  metal- semiconductor  alloys.  As  crude 
approximation,  we  shall  assume  it  to  be  equal  to  that  of  the  pure 
semiconductor.  This  will  permit  us  to  use  the  simplified  potential 
distribution  analysis  of  Section  IV. 

It  would  also  be  interesting  to  simultaneously  investigate 
the  possibility  that  thermoelectric  effects  play  a  part  in  the  forming 
process  (54,  55). 

The  three  most  important  thermoelectric  effects  not  involving 
magnetic  fields  are  the  Seebeck  effect,  the  Thompson  effect,  and  the 
Peltier  effect.  They  are  so  closely  related  that  the  simplest  way 
to  understand  any  one  of  them  is  to  consider  all  three. 

When  a  circuit  composed  of  two  dissimilar  materials  is 
such  that  one  junction  is  at  a  different  temperature  than  the  other, 
(see  Figure  19)  there  will  be  developed,  in  general,  an  E.M.F.U^^- 
The  thermoelectric  power  then  defined  by  Equation  40. 

ai2  =  lim  .  (40) 

6T-0 

AT  is  the  temperature  difference  between  the  two  junctions. 
We  shall  discuss  the  sign  convention  later. 

The  Peltier  effect  may  be  thought  of  as  the  converse  of  the 
Seebeck  effect.  When  current  is  passed  across  the  junction  between 
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two  dissimilar  metals,  the  junction  absorbs  or  liberates  a  quantity 
heat  per  unit  area,  proportional  to  the  current  density  J. 


FIGURE  19  THERMOELECTRIC  EFFECTS. 

In  contrast  to  the  Peltier  and  Seebeck  effects, the  Thompson 

effect  relates  to  a  single  material  rather  than  to  contacting 

materials.  The  Thompson  effect  relates  to  the  rever sable 

absorption  or  liberation  of  heat  which  occurs  when  a  current  flows 

along  a  thermal  gradient.  Equation  4Z  defines  the  Thompson 

dT 

coefficient  T  in  terms  of  the  thermal  gradient  the  current 
density  J  and  the  heat  per  unit  volume  generated. 
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(42) 


At  this  point,  a  word  regarding  signs  is  in  order.  Suppose 
in  Equation  (40)  that  AT  is  always  taken  as  positive.  Then  U^2 
considered  positive  if,  in  an  open  bimetal  circuit  such  as  that 
shown  in  Figure  19,  the  terminal  connected  to  the  high  temperature 
side  is  positive.  Still  referring  to  Figure  19,  2  heat 

withdrawn  from  the  environment  at  terminal  a,  while 
heat  withdrawn  from  the  environment  at  terminal  b.  Thus,  if 
terminal  b  heats  up,  then  ^2^  ^  negative  number;  if  it  cools 

down,  Q2^  is  a  positive  number.  In  Equation  (42),  if  T  is  positive 
and  if  the  conventional  flow  of  current  is  in  the  direction  of 


increasing  temperature,  then  heat  is  absorbed  from  the  environment. 

The  three  thermoelectric  coefficients  are  linked  by 
Equations  (43)  and  (44),  commonly  known  as  the  Thompson 
relations; 


Although  these  relations  were  originally  developed  using 
classical  concepts,  they  have  been  proven  valid  by  quantum 
mechanical  treatments  also. 
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Equation  (44)  permits  us  to  divide  ^  i’^to  two  portions, 
each  of  which  is  dependent  upon  one  material; 


CL  1  ^  ss  0l_  -  CL . 

1,2  2  1 

where 

a  =  r  T  <*’' 

^  .0 

and 


“"l 

T~ 


da 


1 


“  -3T 


(45) 

(46) 


(47) 


It  thus  becomes  meaningful  to  speak  of  the  thermoelectric 
power  of  a  single  material.  Experimentally  this  is  a  great  help, 
since  it  is  often  much  easier  to  measure  thei*moelectric  power 
than  to  measure  the  Thompson  coefficient  directly. 

Unfortunately  for  our  purposes,  neither  theoretical  nor 
experimental  work  has  been  done  on  the  thermoelectric  parameters 
of  germanium  as  heavily  doped  as  that  which  we  are  using.  Johnson 
and  Lark-Horowitz  (56)  considered  the  case  where  the  doping  was 
sufficient  that  ionized  impurity  scattering  was  the  dominant  obstacle 
to  carrier  motion;  however,  they  did  not  consider  degenerate 
material  and  employed  Boltzmann  Statistics.  Their  results, 
together  with  the  experimental  results  of  Middleton  and  Scanlon 
(57),  are  shown  in  Figure  20.  The  knees  of  the  curves  occur  at 
temperatures  where  intrinsic  conduction  becomes  important.  We 
have  seen  that  in  our  case  this  happens  only  at  temperatures 
approaching  the  melting  point  of  germanium. 
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55 

Johnson  derives  an  expression  for  the  Thompson 
coefficient  valid  for  semiconductors  at  all  doping  levels. 
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and,  combining  these  with  (46), 
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(51) 


e  is  the  charge  on  the  electron,  e  is  the  energy,  measured 
from  the  bottom  of  the  valence  band,  fp(6)  is  the  unperturbed 
energy  distribution  function,  f  is  the  mean  free  path,  ^  is  the 
Fermi  level,  and  m  is  the  electron  effective  mass. 

It  might  be  noted  that  is  the  conductivity  of  the  material. 
One  is  tempted  to  use  the  results  of  Section  II  and  assume  it 
constant;  however,  we  must  bear  in  mind  that,  in  Section  11,  we 
were  only  interested  in  the  variation  in  resistivity  compared  to 
the  average  resistivity,  while  here  we  are  interested  in  the 


derivative. 
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It  is  now  necessary  to  obtain  some  of  the  terms  occurring 
in  Equations  (48,  49,  and  50).  The  variation  of  ^  with  temperature 
in  heavily  doped  germanium  was  studied  by  Blakemore  (58);  his 
results  are  shown  in  Figure  21.  Note  that  the  energy  gap  between 
the  bottom  of  the  conduction  band  and  the  top  of  the  valence  band 
narrows  with  increasing  temperature. 

In  order  to  write  the  f^  term,  it  is  first  necessary  to  know 
N(c).  An  expression  which  takes  rough  account  of  the  effect  of  the 
many -valley  E(k)  characteristic  and  of  the  anisotropic  energy 
surfaces  is  given  in  Equation  (52). 

1.T/  \  Z  "1^^^  1/2  v-.o29  1/2 

N(€)  —  ^ L  II  ^  “3.51X10  c 

(52) 

where  m  ||  and  m  j^are  effective  masses  relating  to  the  curvature 
of  the  band  along  two  different  directions. 

The  mean  free  path  I  may  be  calculated  by  means  of  the 
Conwell- Weisskopf  formula  (36). 


I  = 


,222 
k  G  c 
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vNefn 


,^,2  -4 

1+k  €  d  e 
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(53) 


\ 


where  k  is  the  dielectric  constant  of  germanium,  g^  is  the 
permittivity  of  free  space,  g  is  energy,  N  is  the  dopant  density,  d 
is  the  average  distance  between  impurities,  and  e  is  the  electronic 
charge.  The  calculation  was  carried  out  and  the  results  are  shown 
in  Figure  22. 
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*  A/y 

N(e)  =  ^  [2^11  =  3.  51  X  10^*^ 

(52) 

where  m  |  and  mj^are  effective  masses  relating  to  the  curvature 
of  the  band  along  two  different  directions. 

The  mean  free  path  I  may  be  calculated  by  means  of  the 
Conwell-Weisskopf  formula  (36). 
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where  k  is  the  dielectric  constant  of  germanium,  is  the 
permittivity  of  free  space,  e  is  energy,  N  is  the  dopant  density,  d 
is  the  average  distance  between  impurities,  and  e  is  the  electronic 
charge.  The  calculation  was  carried  out  and  the  results  are  shown 
in  Figure  22. 


-69- 


8 


lO 


o 

o 

CM 


o 

o 


o 

o 

o 


8 


<Tt 


Q  5^ 

e 


S 


I 


o 

o 

m 


O 

O 


to 


FIGURE  21  VARIATION.  OF  FERMI  LEVEE  WITH  TEMPERATURE  N  =  ^  v  -in^9  atoms 
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We  can  combine  (52)  with  the  Fermi-Dirac  distribution 


function  to  obtain  f^Ce) 


fo(e)  = 


3.  51  X  10^'^e^^^ 


combining  this  with  (51)  and  simplifying, 
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(55) 


This  was  calculated  using  the  values  of  ^  from  Figure  21 
and  of  f  from  Figure  22.  The  results  are  shown  in  Figure  2  3;  they 
appear  to  be  a  reasonable  extension  of  the  data  in  Figure  20. 

The  slope  of  the  thermoelectric  power  curve,  multiplied 
by  T,  is  the  Thompson  coefficient.  To  a  good  approximation  the 
curvature  of  a(T)  may  be  neglected,  and  we  may  write: 

C  =  -0.16  T  microvolts/°K  (56) 

for  germanium  of  this  doping  level. 

We  shall  now  generalize  Equation  (42)  to  a  three-dimensional 

2 

case  and  add  the  heat  source  so  obtained  to  the  I  R  heat  source  in 
Equation  (34). 
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4V^a 

V.  (k7T)  -  .  .  -  x-°  ° - 5-  +  tT-  VT  =  0 

a^(r+C)(l+^^  ) 


(57) 


The  current  density  may  be  derived  from  the  expressions 
of  Section  IV. 


T  =  -  c.vv  = 


2V  c 
o  s 


(58) 


Next  we  express  V  •  (kVT)  in  terms  of  coordinates 

h^h 

a  I  1. 


=  hAjh,  -k 


Fk  n 


a2(^2^^2)3l^ 


and  finally,  writing  out  the  third  term  in  (57) 
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We  shall  evaluate  the  temperature  profile  along  the  centerline 
(I  =  1).  Then  (57)  becomes 
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V  (i+r) 

This  was  solved  by  nvunerical  computation  for  =1.3,  which 
was  a  value  frequently  employed  in  forming  diodes.  The  thermal 
profile  is  shown  in  Figure  24.  The  liquid- solid  interface  may 


FIGURE  24  TEMPERATURE  PROFILE  FOR  Vo=  1.  3v  .  C-t  =  SOLID  LIQUID  INTERFACE. 
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be  taken  to  be  the  surface  for  which  T(C)  =  700°K,  the  eutectic 
temperature  for  the  germanium,  and  corresponds  to  C  =  1.2. 

It  might  be  noted  that  the  Thompson  heating  is  negligible 

in  comparison  to  the  joule  heating  except  in  the  region  immediately 

adjacent  to  the  interface,  where  the  high  values  of  T  and  dT/dC 

2 

cause  it  to  rise  to  20  per  cent  of  the  I  R  term. 

It  is  often  of  interest  to  be  able  to  calculate  the  thermal 
time  constant  of  a  structure  of  this  sort.  Actually,  there  are 
two  thermal  time  constants  in  a  system  of  this  type.  The  first 
time  constant  is  the  length  of  time  after  the  application  of  a 
voltage  which  would  be  required  to  reach  the  steady  state  temperature 
corresponding  to  that  voltage,  assuming  that  the  initial  rate  of 
change  of  temperature  were  held  constant.  The  second  time  con¬ 
stant  is  the  length  of  time  after  the  removal  of  a  voltage  for  the 
temperature  to  return  to  the  ambient  temperature,  providing  that 
the  rate  of  decrease  of  temperature  remained  equal  to  what  it  was 
immediately  after  the  removal  of  the  voltage.  For  most  practical 
purposes  we  may  regard  the  two  thermal  time  constants  as  equal, 
at  least  to  within  an  order  of  magnitude. 

To  see  why  this  is  so,  let  us  first  agree  to  ignore  the 
Thompson  heating,  an  approximation  which  we  have  seen  should 
cause  an  error  of  20  per  cent  or  less.  Then  we  may  divide  the 
heat  flow  into  two  portions:  that  due  to  the  thermal  gradient  and  that 
due  to  the  ohmic  power  loss.  Let  us  call  the  former  Qrp  and  the 


latter  Qj.  Then: 
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I  “  2  2,^1  ^2, 


and 


=  V.  (kVT) 


(63) 


(64) 


Before  equilibrium  is  established,  the  thermal  continuity- 
equation  must  contain  a  time -dependent  term: 


Qt  +  Qi  -  c 


(65) 


where  C  is  the  specific  heat  of  the  material. 

Suppose  now  that  at  time  t  =  t^  we  apply  a  voltage  pulse 
to  a  "  cold"  diode.  Since,  initially,  T  =  the  ambient  temperature 
throughout,  Q.^,  is  identically  zero,  Q^,  however,  is  independent 
of  temperature  (as  we  showed  in  H)  and  hence  of  time,  provided 
only  that  the  applied  voltage  does  not  vary.  Thus, 

dT 


Qt  +  Qj 
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=  C 


t  =  t 


(66) 


Then  the  thermal  time  constant  relating  to  heating  up, 


which  we  shall  call  t^,  is 


"  9WatT^ 

Here  T£  is  the  final,  steady -state  temperature.  When  a 
steady- state  has  been  reached  Q.p  =  -Qj  and  the  time  derivative 
is  zero.  When  the  voltage  is  removed,  Qj  becomes  identically 
zero  and  Qj  does  not  change  instantaneously,  since  the  temperature 


T£C(/a^)(e^C^)(lH^) 
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does  not  change  instantaneously  at  any  point.  Hence, 
voltage  is  removed  at  t  =  t^, 


Q,j.,  +  Qj  =  Q.J,  =  -Qj 
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and,  defining  the  cooling -down  time  constant  as  t^ 
T.C(Tr^a^)(eV^)(l+^^ 

'2  = - - 

4V  a 
o  3 

C  is  independent  of  temperature  in  the  temperature  range 
of  interest  here,  and  so  is  a  g.  Hence  ^2  ~  Thompson  heating 

were  taken  into  account,  it  would  be  a  factor  in  t^  but  not  in  t^. 

To  get  some  idea  of  the  order  of  magnitude  of  t^,  let  us 

take  the  thermal  time  constant  relating  to  the  solid-liquid 

interface  in  the  case  where  a  =  lO""^  cm,  V  =  1.  3V,  T,  =  700°K, 

o  ’  f  ' 

3 

=  10  ,  ^  =  1,  and  ^  =  l.Z.  C  =  5.47  calories/mole  or 
1. 68  joules/c.  c.  For  these  values  t^  =  8.  3  jisec. 

This  is  a  quantity  which  one  may  check  experimentally. 
Suppose  that  we  have  a  configuration  with  a  particular  value  of  a, 
and  we  apply  a  pulse  of  height  sufficient  to  form  a  good  diode,  but 
with  a  length  which  is  close  to  the  thermal  time  constant.  Then 
we  may  expect  that,  whatever  the  I-V  characteristic  of  the 
resulting  diode  might  be,  a  second  pulse  of  the  same  length  will 
alter  it.  If,  however,  a  pulse  is  applied  which  is  long  compared  to 
the  thermal  time  constant,  the  I-V  characteristic  should  remain 
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relatively  constant  with  succeeding  pulses,  presuming  that  no 
additional  factors  such  as  solid  state  diffusion  enter  the  picture. 

In  the  experiment  a  point  was  lightly  performed  to  a 
freshly  cleaved  germanium  wafer  in  such  a  manner  that 

_  3 

a  »  10  cm  (the  technique  for  accomplishing  this  will  be 
discussed  in  Section  VI).  The  current  at  Vp  =  400  mv  was 
taken  as  a  measure  of  the  I-V  characteristic.  This  quantity  is 
plotted  against  the  ntimber  of  pulses  for  various  pulse  len.gths 
in  Figure  25.  The  wire  was  torn  loose  from  the  semiconductor 
before  each  run.  Note  that,  for  a  pulse  length  of  6.  5  (jisec,  the 
value  of  I^QQ  changes  considerably  from  one  pulse  to  the  next, 
whereas  for  the  pulses  long  compared  to  the  predicted  time 
constant,  the  variation  was  slight. 
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VI.  FACTORS  INFLUENCING  JUNCTION  AREA 

There  are  three  factors  effecting  the  junction  area:  the 
pressure  on  the  wire,  the  geometry  of  the  tip  of  the  wire,  and  the 
pulse  supply  impedance. 

It  is,  of  course,  impossible  to  touch  the  wire  to  the 
semiconductor  and  still  have  zero  force  upon  it.  Assuming  then 
that  the  wire  is  bearing  down  on  the  semiconductor  with  some 
nonzero  force,  and  that  a  pulse  of  sufficient  length  from  a  source 
of  nonzero  impedance  is  applied,  the  tip  of  the  wire  melts,  spreads 
out,  and  permits  more  current  to  flow.  Assuming  that  the  pulse 
is  of  siafficient  duration,  this  continues  until  one  of  two  things 
happens:  either  the  IR  drop  in  the  pulse  supply  causes  the  applied 
voltage  to  fall  below  the  critical  value,  or  enough  of  the  wire 
melts  to  release  the  downward  pressure.  In  the  former  case,  we 
are  left  with  the  aluminum  or  other  point  material  spread  out  over 
the  surface,  while  in  the  latter  the  voltage  is  high  enough  to  cause 
the  formation  of  a  rather  extensive  alloy  region  beneath  the  point 
extending  into  the  body  of  the  semiconductor. 

We  shall  see  in  the  next  section  that  the  latter  situation 
may  be  expected  to  lead  to  good  tunneling  junctions.  It  follows 
from  this  that,  even  with  supplies  of  moderately  high  impedance, 
we  should  be  able  to  obtain  txinnel  diodes,  providing  that  the  area 
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of  the  contact  is  small  enough.  It  also  follows  that,  for  a  given 
precontact  point  geometry,  the  allowable  source  impedance 
goes  up  as  the  contact  pressure  goes  down. 

This  was  checked  experimentally  and  found  to  be  true. 

Wires  were  brought  in  to  contact  with  the  crystal  at  pressures  and 
with  point  geometries  which  normally  required  0.  1  ohm  pulse 
sources  to  yield  tunnel  diodes.  It  was  tacked  by  means  of  a  100- 
1000  ohm  source,  with  pulse  lengths  sufficient  to  cause  the 
initial  current  level  of  the  (non-tunnel)  diode  to  reach  some 
arbitrary  level.  The  portion  of  the  fixture  holding  the  wire  was 
then  moved  backwards  some  amount  determined  by  trial  and  error 
and  dependent  upon  the  amount  of  springiness  in  the  wire:  usually 
between  0.  0001'*  and  0.  OOl’*.  A  pulse  from  a  source  having  an 
impedance  of  between  1  and  100  ohms  was  applied,  and  usable 
tunnel  diodes  of  relatively  low  peak  currents  but  good  peak-to- 
valley  ratios  were  obtained.  (Figure  26).  Without  this  backing- 
off  process,  which  of  course  reduced  the  pressure  on  the  contact, 
pulses  of  this  impedance  level  rarely  formed  good  tunnel  diodes. 
Needless  to  say,  if  the  wire  is  backed  off  too  much,  the  forming 
pulse  simply  causes  the  diode  to  open,  or  else  the  very  delicate 
tacked-on  contact  breaks  before  pulsing.  The  permissible 
impedance  of  the  pulse  supply  may  also  be  raised  by  taking  extreme 
care  to  make  the  irdtial  contact  as  light  as  possible. 
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(a ) 


(b) 


FIGURE  26  EFFECT  OF  HEAVY  PRESSURE  (A)  AND  LIGHT 
PRESSURE  (B)  IN  THE  FORMING  QF  DIODES 
WITH  THE  SAME  INITIAL  JUNCTION  AREA.  IN 
BOTH  CASES,  FORMING  PULSE  WAS  3,0  VOLTS, 
10  OHMS,  AND  130  MICROSECONDS. 
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vn.  THE  IMPURITY  PROFILE  OF  THE  FINISHED  DIODE 

The  theoretical  and  experimental  results  up  to  this  point  have 
suggested  two  main  types  of  diode  structures.  The  first  type,  illus¬ 
trated  in  Figure  14,  is  characteriited  by  a  metal- semiconductor  sur¬ 
face.  This  type,  typically  resulting  from  high-impedance  forming 
network  as  described  in  Section  VH,  will  henceforth  be  referred  to  as 
planar  configuration.  The  second  type,  which  we  shall  term  a  spher¬ 
oidal  type,  is  exemplified  bythe  diode  depicted  in  Figure  18  and 
Figure  24.  It  is  characterized  by  a  region  of  metal- semiconductor 
eutectic  extending  down  into  the  body  of  the  semiconductor  and  up¬ 
wards  into  the  wire. 

This  latter  picture  of  the  diode  is  substantiated  by  examination 
of  the  regions  left  when  the  wire  was  pulled  away  after  forming. 

(Figure  11).  Ignoring  the  shear  area  on  the  periphery  (presumably 
due  to  thermal  strains  set  up  in  the  germanium  lattice),  the  area  where 
the  wire  was  pulled  away  appears  to  be  a  smooth  spheroidal  surface. 
Estimates  of  its  proportions  in  the  x-direction  were  possible  because 
of  the  narrow  depth  of  field  of  the  metallographic  objective.  By  cal¬ 
ibrating  the  amount  of  motion  needed  to  focus  first  at  the  base  of  the 
hemispheroid,  then  on  its  top,  a  z -dimension  extending  above  the 
surface  of  the  semiconductor  somewhat  less  than  half  the  radius  of  the 
contact  was  found.  The  base  of  the  spheroid  was  somewhat  less  than 
this  distance  below  the  surface.  This  is  roughly  what  we  would  expect 
if  the  surface  were  the  eutectic -aluminum  interface.  Because  of  the 
differential  expansion  upon  cooling  between  eutectic  and  aluminum,  and 
possible  also  because  of  a  difference  in  lattice  configuration,  this 
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interface  is  very  likely  a  stress  concentration  plane,  and  might 
very  well  be  the  first  place  to  fail  when  the  wire  is  pulled  away. 

It  has  been  the  experience  of  the  Bell  Laboratories  group 
(2,13.14),  as  well  as  our  own  experience,  that  low-impedance  pulse 
supplies  (which  we  have  seen  lead  to  spheroidal  configurations)  are 
necessary  to  form  tunnel  diodes.  We  are  now  in  a  position  to  explain, 
at  least  qualitatively,  why  this  should  be  so. 

It  was  mentioned  earlier  that  diodes  of  this  sort  depended  upon 
a  process  similar  to  dot-alloy  diodes  for  junction  formation.  Let  us 
now  examine  this  notion  a  little  more  closely. 

Alloy  junctions  (59)  are  based  upon  the  principle  that,  when  a 
melt  consisting  of  a  semiconductor  and  a  dopant  metal  cools  in  contact 
with  germanium,  it  is  possible  for  the  germanium  to  regrow  expitaxially 
on  the  solid  germaniiim,  incorporating  within  it  small  qviantities  of  the 
dopant  metal.  The  alloy -semiconductor  combination  in  a  dot-alloy  unit 
cools  by  two  mechanisms:  conduction  through  to  semiconductor  to  a  heat 
sink,  and  convection  and  radiation  from  the  metal.  The  former  ordinarily 
dominates  the  latter,  although  sometimes  auxiliary  heaters  abovfe  the 
wafer  are  recommended  to  further  slow  down  the  cooling  of  the  top  of  the 
melt.  If  the  rate  of  cooling  of  the  external  part  of  the  melt  begins  to 
become  comparable  to  that  of  the  allow- semiconductor  interface,  a 
situation  can  result  in  which  nucleation  of  the  regrowing  germanivun  takes 
place  first  at  point  jutting  into  the  alloy.  This  leads  to  a  jagged 
junction  and  often  poor  diode  characteristics. 

With  a  point  contact  diode,  the  situation  is  even  worse  with  res¬ 
pect  to  the  formation  of  high-quality  p-n  jxmctions.  Here  the  portion  of 
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the  melt  away  from  the  germanivim  is  coimected  to  the  point,  which 
is  a  heat  sink  having  a  thermal  conductivity  many  times  that  of  the 
semiconductor.  Thus,  in  a  planar  type  diode,  the  thin  molten 
region  faces  two  parallel  heat  paths,  the  one  corresponding  to  the 
point  having  a  far  higher  thermal  conductance  than  that  correspond¬ 
ing  to  the  semiconductor.  Under  these  conditions  it  is  to  be  expec¬ 
ted  that  the  nucleation  will  take  place  from  the  outside  in,  with  little  or 
no  nucleation  onto  the  solid  germanium,  the  germanium  rejected  by  the 
cooling  alloy  being  precipitated  in  isolated  inclusions.  Under  these 
conditions  the  re  suiting  diode  will  be  a  poor -quality  metal- semiconductor 
type.  This  has  indeed  been  the  case  whenever  diodes  have  been  pulsed 
from  the  high -impedance  sources  which  we  have  seen  encourage  the 
planar  configuration. 

The  situation  which  prevails  in  the  case  of  apronounced  spher¬ 
oidal  configuration  such  as  that  found  in  Section  Vis  considerably  more 
favorable  to  junction  formation.  Here  the  thermal  resistance  which 
the  alloy -semiconductor  interface  sees  is  that  of  the  alloy  between  the 
interface  and  the  wire. 

To  see  the  importance  of  this  effect  let  us  calculate  the 
thermal  resistance  looking  in  both  directions  from  the  interface. 

For  our  purposes,  we  may  achieve  svifficient  accuracy  by  asstaming 
the  thermal  conductivities  of  the  semiconductor  and  the  alloy  are 
independent  of  temperature  and  equal  to  K  and  K  ,  respectively. 

S  a 

Let  =  the  value  of  ^  corresponding  to  the  interface,  and  At  = 

X  s 

the  temperature  difference  between  the  junction  and  the  ambient 
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(T(oo)).  Then,  solving  Equation  (24)  >vith  T  substituted  for  V 

and  K  substituted  for  a  ,  we  get 
s  s 

A  T 

T  - 2 — ^  4J+  T  ambient  .  (70) 

1  — tan"  Ct 
TT  I 

The  heat  flow  away  from  the  junction  will  be 


H  =  r  K  V  T  da= 

BB  B 

area  of 
junction 


2AT  K 
_ s  a 

iTa(l  -^an"^^^) 


/ 


da 


area 
junction 


ofV^ 


(71) 


The  germanium  aluminum  eutectic  is  70  per  cent  alximinum, 
80  we  would  expect  the  alloy  to  project  above  the  plane  of  the 
semiconductor.  However,  as  a  conservative  estimate  let  us 
assume  that  the  alloy -aluminum  interface  lies  at  ^  =  0.  Then, 
assuming  that  we  have  a  planar  heat  source  of  temperature 
6T  +  T  ambient  at  ^  and  that  the  alloy -aluminum  interface 
is  at  the  ambient  temperature, 

AT 

T  =  T  ambient  +  - ? —  tan  ^  .  (72) 

Tsin\ 

Then  the  heat  flow  away  from  the  junction  in  the  direction 
of  the  aluminum  is  then 


H  =  r  K 

as  J  a 

area  of 
junction 


VT  da  = 


AT  K  „ 

8  a  r 

aTan 

area  of 
jvinction 


da 

(e^+Ci)(i+Cj) 


and  the  ratio  of  the  thermal  resistivities  is 


(73) 
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2K  tan'^T 

S  X 


.  (74) 


iraCl  tan'^t,  ) 
TT  X 


In  a  planar  structure  0  and  the  ratio  of  resistances 
becomes  infinite.  Hence,  if  we  take  an  electrical  analogue  and 
picture  the  layer  of  molten  eutectic  at  the  alloy- semiconductor 
interface  as  a  charged  capacitor,  most  of  the  heat  will  flow  towards 
the  wire,  and  we  would  expect  the  portion  nearest  the  wire  to  cool 
down  first. 

On  the  other  hand,  if  we  cite  our  example  of  part  V  for 
a  heavily  formed  diode,  ®  ^  and: 


K 

a 

K“ 


(75) 


a  s 

The  alloy,  being  a  highly  disordered  structure,  will  ordinarily  have 
a  considerably  lower  thermal  conductivity  than  the  semiconductor, 
hence  R  <  R  ,  an  appreciable  quantity  of  heat  flows  through  the 
semiconductor  and  the  heat  of  fusion  of  a  layer  of  eutectic  will  have 
had  time  to  have  been  removed  from  a  layer  bounding  the 
semiconductor  before  the  freezing  front  from  the  aluminum  side 
has  traversed  the  melt,  and  good  regrowth  and  usable  tunnel  diodes 


are  the  result. 
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The  objection  might  be  raised  that  initial  thermal  gradients 
were  ignored  in  this  analysis.  However,  this  also  was  a  conservative 
approximation,  since  initially  the  alloy- semiconductor  interface 
looks  at  a  steep  thermal  gradient  in  the  semiconductor  direction 
and  no  gradient  at  all  in  the  alloy  direction  (the  whole  melt  is 
presumed  to  be  at  the  eutectic  temperature). 

At  this  stage  it  might  be  well  to  mention  very  briefly  solid 
state  diffusion  of  the  point  material  into  the  unmelted  crystal. 

The  diffusion  constant  D  which  relates  the  flux  of  dopant 
atoms  to  the  concentration  gradient  is  characterized  by  a  multi¬ 
plicative  constant  and  an  activation  energy  AH. 

D  =  exp  [ah/Rt]  (76) 

where  R  is  the  gas  constant.  For  aluminvim  in  germanium,  there 
f;ippears  to  be  no  data  on  AH  on  D^.  To  get  an  idea  of  the  order 
pf. magnitude  of  the  effect  at  the  temperature  which  we  used  in  the 
example,  we  shall  average  the  values  of  the  neighboring  column  HI 
elements  in  germanium  (one  should,  however,  have  no  illusions 
regarding  the  validity  of  such  a  procedure).  This  yields 

•j 

AH  «  80  kilocalories  and  D  =  1000  cm'^/sec.  At  T  =  700°K, 

o  ' 

-18 

this  means  D  is  on  the  order  of  magnitude  of  10"  .  If  we  ignore 

the  curvature  of  the  interface,  we  get  for  the  gradient  distribution 

N(+,  z)  =  Nt  erfc  — - -  (77) 

2  VD+1 
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where  N  is  dopant  density,  Nj  is  the  density  at  the  interface 

(assvuned  constant)  and  t  is  the  time  from  the  application  of  heat. 

20 

The  solid  solubility  of  aluminum  in  germanium  is  about  10 

atoms/cc.  If  we  are  interested  in  the  period  necessary  for  an 

18 

appreciable  number  (say  10  /cc)  of  dopant  atoms  to  diffuse  100 

Q 

Angstroms,  then  we  derive  a  value  of  t  equal  to  10  seconds. 
With  all  due  account  of  the  inadequacy  of  our  estimates  of  AH 
and  D^,  it  seems  unlikely  that  pulsing  formats  of  less  than  a 
second  will  give  rise  to  solid-state  diffusion  effects  at  these 
temperatures. 
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Vin.  CONCLUSIONS 

Before  stating  the  conclusions  reached,  it  might  be  well 
to  mention  some  factors  which  were  not  considered  in  this  study. 

First,  we  considered  a  range  of  dopants  and  temperatures 
where  intrinsic  conduction  was  never  important.  However,  it  is 
quite  possible  that  this  assumption  would  be  invalid  in  the  case  of 
backward  diodes  and  some  tunnel  diodes. 

It  was  always  assumed  that  electrons,  when  they  were 
scattered,  had  their  momentum  completely  randomized.  However, 
it  is  possible  that  with  the  extremely  high  fields  which  exist  in  the 
neighborhood  of  the  point,  **hot  electrons"  might  play  some  role 
in  transport  processes. 

The  motion  of  ionized  impurity  atoms  under  the  influence 
of  high  electric  fields  and  high  temperatures  has  received  some 
attention  In  recent  years  and  might  play  some  role  in  determining 
the  impurity  profile.  Unfortunately,  certain  experimental  advan¬ 
tages  have  caused  most  of  the  work  in  this  area  to  be  carried  out 
with  lithixim  as  the  diffusant,  which  makes  quantitative  work  with 
ni-V  dopants  difficult. 

Peltier  heating,  which  received  little  more  than  a  passing 
mention  here  since  it  took  place  on  the  alloy  side  of  the  junction, 
might  be  important  in  diodes  which  are  reversed  pulsed. 


Interaction  between  dopants  was  ignored;  it  was  assumed 
that  the  diffusion  of  aluminvim  into  the  arsenic -doped  material  dom 
inated  the  converse  process  because  of  the  very  much  higher 
concentration  of  aluminum.  For  an  argument  which  indicates 
that  this  might  not  be  the  case,  the  reader  is  referred  to  the 
article  by  McCaldin  (61). 

En\jmerated  below  are  the  contributions  which  the  author 
feels  this  work  has  made  to  the  literature  of  electrical  engineering 

1)  It  was  shown  that  the  rather  simple  model  used  by 

Sim  (1)  and  Torrey  and  Whitmer  (3)  was  adequate  to  give  accurate 
numerical  estimates  of  the  minimum  voltage  necessary  to  give 
any  forming  effect,  provided  that  one  takes  into  account  the 
temperature  variation  of  thermal  conductivity  and  the  current 
limiting  effect  of  the  metal- semiconductor  junction.  This  is 
the  first  aoad  only  study  of  forming  which  takes  these  two  factors 
into  account. 

2)  The  thermoelectric  power,  Thompson  coefficient, 
and  resistivity  of  extremely  degenerate  germanium  at  elevated 
temperatures  was  calculated.  Although  the  computation  was 
straightforward  and  involved  well-known  techniques,  the 
parameters  for  this  degree  of  doping  and  range  of  temperatures 
have  not  previously  been  reported  in  the  literature. 
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3)  It  is  shown  that  thermoelectric  effects  can  indeed  be  a 
significant  factor  in  forming. 

4)  It  is  shown  that  a  simple  planar  model  of  the  sort  used 
by  previous  investigators  does  not  lead  to  a  situation  in  which 
single -crystal  regrowth  {which  a  considerable  number  of 
thermoelectric  probing  experiments  have  indicated  occurs)  is 
possible.  A  model  of  the  forming  process  in  which  a  alloy  region 
extends  into  the  body  of  the  semiconductor  is  developed,  and  fits 
well  with  relevant  experimental  data.  In  particular,  the  thermal 
time  constant  predicted  by  the  theory  appears  to  be  essentially 
correct,  and  the  configuration  of  the  contact  after  the  point  is 
torn  away,  which  might  otherwise  be  somewhat  puzzling,  fits 
directly  into  the  model.  Furthermore,  it  is  shown  that  this 
picture  leads  to  a  heat  flow  pattern  in  which  conditions  for 
satisfactory  regrowth  are  met.  The  prediction  of  this  theory  , 
namely  that  a  low  impedance  pulse  supply  is  necessary  for  good 
diode  forming,  agrees  with  our  own  observations  and  those  of 
other  workers  in  the  field, 

6)  Fitially,  it  is  shown  that,  at  least  under  a  rather  typical 
set  of  conditions,  solid  state  diffusion  may  be  ignored  entirely, 
unless  forming  times  are  exceedingly  long.  This  also  agrees 
with  empirical  observation. 
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